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We investigate cosmic-ray antiprotons emitted from the galactic primordial black holes in the 
Randall-Sundrum type-2 braneworld. The recent results of the BESS antiproton observation implies 
the existence of exotic primary sub-GeV antiprotons, one of whose most probable origin is PBHs in 
Our Galaxy. We show that the magnitude of antiproton flux from PBHs in the RS braneworld is 
proportional to negative power of the AdS radius, and immediately flnd that a large extra-dimension 
can relax upper-limits on the abundance of the galactic PBHs. If actually there are more PBHs than 
the known upper-limit obtained in the pure 4D case, they set a lower bound on the size of the extra 
dimension above at least 10^" times 4D Planck-length to avoid inconsistency. On completion of the 
numerical studies, we show that these constraints on the AdS radius is comparable to those obtained 
from the diffuse photon background by some of the authors in the previous paper. Moreover, in 
the low accretion-rate case, only antiprotons can constrain the braneworld. We show that we will 
detect signatures of the braneworld as a difference between the flux of the antiprotons predicted in 
4D and 5D by future observations in sub-GeV region with a few percent precision. 

PACS numbers; 98.80.Cq, 04.50.-|-h, 04.70.Dy, 96.40.-z 



I. INTRODUCTION 



The primordial black holes (PBHs), which are thought 
to be formed from density fluctuations during the 
radiation-dominated era of the early Universe P, Hi 
HilEIEQ (for other mechanisms, see also Q), have long 
been a useful tool to investigate the early Universe it- 
self owing to their radiative nature discovered by Hawk- 
ing [i[l3[ll|ll|lll3E|ll|l3- A cosmological 
model under investigation in this paper is the so-called 
RS2 braneworld suggested by Randall and Sundrum. In 
their scenario, our Universe is a 4D Lorentz metric hy- 
persurface called brane, embedded in a Z2 symmetric 5D 
anti-de Sitter (AdS) spacetime called bulk [Tg- RS2 as- 
sumes that the number of branes is only one and its ten- 
sion is positive. The bulk has a typical length scale that 
originates from AdS curvature radius, I, which draws a 
short distance boundary on the brane; 4D Newton law is 
recovered as far as we consider a sufficiently larger scale 
than I, while deviation from ordinary gravity appears 
below it [11 Ell. The present table-top upper bound 
on the bulk AdS radius I was obtained by measuring 
short-distance gravitational force as I < 0.2 mm poll2l| . 
From such a point of view, PBHs again retain their sta- 
tus. Actually, in the framework of the braneworld, it 



has been revealed that the nature of black holes with 
small masses highly differs from the ordinary 4D one. 
In particular small black holes with Schwarzschild ra- 
dius rs ^ I behave as essentially 5D objects and their 
radii arc significantly stretched compared with 4D BHs 
I22, 23]. It means that the temperature of Hawking ra- 
diation drops down. We naturally expect that even such 
small black holes were produced by gravitational collapse 
in the early universe |2J, [23 (there have been discussions 
on the black hole formation in respect of the particle 
scattering processes even in the early Universe p3. l2^ ). 
Thus we have to revisit the particle production by pri- 
mordial black holes taking into account how non-4D char- 
acter modifies it. In addition, the RS2 scenario implies 
a nonstandard high-ener gy r egime before the standard 
radiation-dominated era j27l I28L I29I l30l | (see, e.g., 
and references therein) . It was recently argued by Gue- 
dens, Clancy, and Liddle and Majumdar j2J,|33,|33| that 
through the brane high-energy phase accretion of sur- 
rounding radiation onto PBHs is so efficient that they 
grow in their masses on the contrary to the standard 
scenario. Actually it was already revealed that the high- 
energy diffuse photon background is highly modified in 
the braneworld [34l ISSj . and it can in turn constrains 
braneworld according to the amount of PBHs J5J . Even 
in other various contexts, braneworld PBHs are gradually 
extending their status |36L l37l| . 
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We extend above-mentioned strategy in this paper by 
incorporating cosmic-ray antiproton observations. Be- 
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cause PBHs are thought to be emitting both particles and 
antiparticles without any distinction, we can search for 
dark matter PBHs with the aid of antimatter detectors. 
In fact, one of the most precisely measured is antiproton 
among various antimatter species in cosmic rays. High- 
precision measurements of antiprotons around GeV have 
been successively performed since mid-1990s by several 
collaborations such as Balloon-borne Experiment with a 
Superconducting Spectrometer (BESS) HI US El , 
Isotope Matter Antimatter eXperiment (IMAX) |42j| . 
Cosmic AntiParticlc Ring Imaging Cerenkov Experiment 
(CAPRICE) and so on. There are even observations 
in high energy regions over 3 GeV |3 US ■ It is known 
that antiprotons can set nearly the same constraint on 
the PBH abundance as the diffuse gamma-ray in the 
4D case ^3 QS- Therefore in this paper, we ex- 
amine antiproton emission from braneworld PBHs and 
see the spectral shape of the antiproton flux. Combined 
with the observational data, they give quantitative con- 
straints on the amount of PBHs and, in the context of 
the braneworld, in turn set bounds on the AdS curvature 
radius I. 

For propagation of cosmic-ray antiprotons in the 
Galaxy, we take a recent version of the diffusion model 
lilllflllllEllElllllllIll- Unfortunately, the effect 
of the solar wind is crucial for the sub-GeV antiproton 
flux. Although such a effect, called solar modulation, has 
been under successive measurement, there remains a dif- 
ficulty in incorporating it into the model. Hence we wish 
to use data as little affected by the solar activity as pos- 
sible^herefore choose BESS (slliSlil and CAPRICE 
|43l 1451 , which were taken during the recent solar mini- 
mum. 

The contents of this paper are as follows. First we 
briefly review in Sec. ^ the cosmology and behavior 
of evaporating black holes in the RS2 braneworld. In 
Sec, mil formulations to describe PBHs in the brane Uni- 
verse are shown. We analyze the dependence of the an- 
tiproton flux on the braneworld parameters in Sec. IIVI 
In this section, we also show numerical results for a wide 
range of braneworld parameters to give quantitative con- 
straints on the PBH amount and themselves. Discussions 
and conclusions arc summarized in Sec. 



II. COSMOLOGY AND BLACK HOLES IN THE 
BRANEWORLD 

A. Braneworld cosmology 

Hereafter we basically use 4D natural units h = c = 
fee = 1 throughout this paper. Newton's gravitational 
constant is also taken to be unity and we define 4D Planck 
mass as M4 = G^^/^ = 1. In these units, other 4D 
Planck-scale quantities, defined as I4, ti, and T4, are also 
all unity, whereas they will usually appear in expressions 
to clarify their differences from higher dimensional funda- 
mental scales. A subscript "5" expresses five-dimensional 



quantity, such as 5D Planck mass M5, and so on. 
We start with the fivc-dimensional Einstein equation 
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(5) 



(2.1) 



where the energy-momentum tensor is expressed in the 
Gaussian normal coordinates as [l9j | 



T — 



'A55p^ + S^^S{y), Sf_,^ = Xq^u + T^u, (2.2) 



with the metric on the brane q^^- Cosmological expan- 
sion on the brane is governed by the reduced 4D Einstein 
equation |T^ . 
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(2.3) 



Mt = -^, (2.4) 
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(2.5) 

The standard perfect-fluid energy-momentum tensor on 
the brane t^^ leads to the effective 4D Friedmann equa- 
tion. In this framework, the expansion-law is modi- 
fied from the pure 4D case due to the existence of the 
typical energy scale on the brane, the brane tension A 
EllSllSllillsS. The result is 



Pkk 



A4 
3 



K 
7^' 



(2.6) 



where a is the scale factor, H is the Hubble constant, 
p and pkk are the energy density of ordinary matter 
field and of the dark radiation, respectively, and K is 
the intrinsic spatial curvature of the brane. The energy 
conservation law is the same form as the 4D one: 



p+'iH{p + p) ^ 0. 



(2.7) 



For complete derivation of the effective Friedmann equa- 
tions, see, e.g., jsj and references therein. 

As for the dark radiation term, it was argued that this 
quantity must be sufficiently smaller than ordinary radia- 
tion at the time of nucleosynthesis 30, 54]. Thus we can 
simply drop this term. Additionally, based on WMAP 
observation j53|, we can set K = Q and require A4 to 
be extremely small. Using the de finition of the AdS cur- 
vature radius I = y^5M|/|4^rA5| , we obtain the relation 
between the geometrical scale and the brane tension: 



I 



3_Ml 



(2.8) 



Therefore a following simple relation between the five- 
dimensional and four-dimensional Planck mass holds: 



Ml 



IMl 



(2.9) 
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As a result, the Friedmann Eq. (|2.6() becomes 



cosmological quantities in the early Universe for later use: 



This expression differs from the ordinary 4D one with 
respect to the existence of a p-square term. This term is 
only relevant when the radiation energy is sufficiently 
denser than the brane tension, namely in the earliest 
era of the Universe. Then we immediately understand 
that the cosmological evolution in the early Universe is 
split into two phases; the earlier one is the nonstandard 
high-energy regime in which evolution of the scale fac- 
tor a{t) obeys a{t) cx i^/"*, and the latter is the standard 
radiation-dominated Universe, with a{t) cx t^/^. The for- 
mer is often called the p-square phase. 

In order to obtain whole profile of the scale factor, 
WMAP data, such as io = 13.7 Gyr, h = 0.71, n^fih"^ = 
0.135, r^A.o — 0.73, and Zeq = 3233, are used. The time of 
matter-radiation equality teq is calculated to be 72.6 kyr 
[35l |. There is the analytic solution through the whole 
radiation-dominated era: 



a{t) = Oe 



il/4(i + iJl/4 



1/2 



(2.11) 



where we denoted the time when transition occurs as = 
1/2. The asymptotic forms of the profile is determined 
as 



= Rait) ^ it, M^{t)^ Mil- 



Pit) 



for t < ic, 
3Af2 



Rn{t) = 2t, Mn{t) = Mjt 



for t > tc 



B. Five-dimensional black holes 



(2.17) 



In the braneworld scenario, a massive object localized 
on the brane within typical length scale much shorter 
than ^ I should see the background space-time as an 
effectively flat 5D Minkowski space .59,| (for possible 
deviations and other exotic features, see |60l. IgU Ib^. IgM 
0,113). Therefore a sufficiently small black hole is no 
longer ordinary 4D one, but is approximately described 
as 5D solution^. If that black hole has no charge and 
is not rotating, the solution must be 5D Schwarzschild. 
In terms of asymptotically flat spherical coordinates, the 
5D Schwarzschild metric is written down |66l | and its 4D 
projection is 



1 



C{r) 



(2.18) 



with 



a(i) = «ea — for t<^tc 



-'eq 1/2 

'«lTl72 



for tc^it < teq. 



(2.12a) 
(2.12b) 



eq 



In terms of the Hubble constant 

2t + tc 



4t(i + tc) ' 



(2.13) 



other cosmological quantities, energy density p(t), Hub- 
ble radius i?H(i), and horizon mass MH(i), are expressed 
as 



Rait) = 



//-2 r 



3Al^ 
16^ 
1 



^AH{t)Hl + 1-1 



AH{t)Hl 



^m{t)Hl + 1-1 



(2.14) 
(2.15) 
(2.16) 



Hereafter we denote the horizon mass at the transition as 
Mc = MH(tc) = (16/27)Z ~ I, and call it transition mass. 
We explicitly write down the asymptotic expressions for 



CM^i-(Il)' 



(2.19) 



where rs is the 5D Schwarzschild radius and dD,2 is the 
line element on unit 2-sphere. In our framework the 
Schwarzschild radius in the five-dimensional spacetime 
is related to the black hole mass M and five-dimensional 
fundamental mass scale M5 as [i^ 



8 M 



1/2 




(2.20) 



while in the ordinary 4D case rs — 2{M /M4)Ia- On the 
other hand, the Hawking temperature of a Schwarzschild 
black hole is determined by space-time periodicity in the 
direction of imaginary time. For a higher-dimensional 
Schwarzschild black hole, the temperature is given in 



^ The emergence of this topological change is not specific for the 
Randall and Sundrum's braneworld scenario, but a generic fea- 
ture for higher-dimensional theories that have some compactifi- 
cation length scales in a broad sense. We will briefly discuss this 
issue in the last of this paper. See Sec. \V\ 
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terms of its radius as Tr = (D - 3)/47rrs |2a,|6g. Thus 
the five-dimensional Hawking temperature is 



Th 



f 



27r7'; 



3 
32^ 



-1/2 



M 



'1/2 



jr) T„ (2.21) 



while inthe4D caseTn = l/47rrs = {II^'k){M /A-Uy'^Ti 



III. FORMULATIONS 



Formation 



It has been argued that in the early Universe primor- 
dial black holes are formed due to gravitational collapse 
caused by density perturbation 0,12, H, IS- We here 
simplify the situation; when a superhorizon-scale density 
perturbation with adequate amplitude enters the Hubble 
horizon at t = ti, it instantaneously collapses into a black 
hole with some fraction of the horizon mass. Hence the 
initial mass of a PBH denoted Mj is given by 



Mi = /MH(ii), 



(3.1) 



where the parameter / expresses above mentioned 
fraction 0, |63 • Utilizing the explicit forms of horizon 
mass in asymptotic regions, i.e., Eq. H2.17|l . we can write 
down the relations between the time of formation and 
initial mass for two types of PBH as 



Mi = i&frhf 

= A 



for ti <C tc 
for ti ^ tc 



(3.2a) 
(3.2b) 



Note that Mi < /Mc corresponds with 5D BHs while 
Afi > fMc with ordinary 4D ones. 



B. Accretion 

Matter field localized on the brane sees a black hole as 
an approximately disk-shaped object with the effective 
Schwarzschild radius [681 



l/(D-3) 



£»-3 



1/2 



(3.3) 



For D — b and £> = 4, rs,off — 2rg and rs^cff = 
(3-\/3/2)rs, respectively. 

Consider an infinitesimal time interval [i,t-|-dt]. Then 
the gain of black hole mass due to absorbing surrounding 
radiation within dt is described as 113, 1331 



dM\ 
"dT j. 



(3.4) 



where p{t) is energy density of radiation. Here we in- 
troduced a parameter F , which determines the accretion 
efficiency] if the geometrical optics approximation holds 



well, namely an incidental relativistic particle can be re- 
garded as a collisionless point particle and its spin can 
be ignored, F should be almost unity, otherwise smaller. 
Due to the lack of knowledge of this issue, we are forced 
to treat F as a free parameter throughout discussions. 

A significant point is here. In the 4D case, the above 
differential equation does not have a solution which de- 
scribes increasing mass. However, in the branc high- 
energy era, it does give a growing solution; accretion oc- 
curs in the braneworld scenario. Substituting relevant ex- 
pressions for the brane high-energy phase into the above 
differential equation, the resultant form becomes 



/dM 



2F M 



for t < tr 



(3.5) 



Integrating it gives the mass evolution formula |32, |3 



M{t) 



2F/iT 



Mi for t<tc. 



(3.6) 



We define the primordial mass Mp of 5D PBHs at the 
end of accretion as Mp = M{tc). Readers should not 
confuse with Mp. For 4D PBHs, = M-,. 



C. Evaporation 

The differential particle emission rate from a D- 
dimensional static uncharged black hole with mass M 
is given by Hawking's formula: 



dt 



D-l 



9j 



(Tj{M,E) dk 
exp{E/TH) ± 1 (27r)^- 



(3.7) 



where subscript j indicates particle species, gj the num- 
ber of internal degrees of freedom of the particle, E = 

^k'^ + the energy including rest mass, aj the absorp- 
tion cross section, the sign in the denominator expresses 
which statistics the particle obeys, and 7h is the Hawking 
temperature. 

The absorption cross section aj is an essential part of 
the famous quantity named greybody factor. It deter- 
mines how far the radiation spectrum deviates from that 
of blackbody. Although black holes under our consider- 
ations are higher-dimensional, ordinary particles will be 
emitted by them into the four-dimensional brane, and 
once the Hawking temperature of a black hole is fixed, 
the emission processes are almost unchanged from the 
4D case. Therefore we can adopt the known energy and 
spin-dependence of absorption cross section aj obtained 
in the 4D case |l^ . It moderately oscillates in the 
lower energy region, but eventually approaches the geo- 
metrical cross section Trrg for any particle species in 
the high-energy regime. 

The life time of a black hole is calculated from the 
summation of all particles' emission rate. Using the for- 
mula H3.7|l under relativistic limit, the total mass decreas- 
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ing rate of a D-dimensional Schwarzschild black hole be- Note that for small five-dimensional black holes we must 

take AcSA = '^''s.off.s- fmattor is the degree of free- 



comes 



cvap 



'dEdf 



-dE 



3 



3 

oo 

1 



D~l 



expiE/Ta) ± 1 (27r)^-i 

A,s,D^D-3/[iD - 2)nD-2] 



dE 



(2^) 



exp{E/Tii) ± 1 
d^;. (3.8) 



Here we assumed that the absorption cross sections aj 's 
are independent of the kind of particle species and en- 
ergy E under the limit. Integration gives the Stefan- 
Boltzmann's law of evaporating black holes: 



/dM 

\~dr 



(3.9) 



cvap 



where we introduced some quantities as follows, go is an 
effective degree of freedom given by contributions from 
bosons and fermions: 



■)D-1 



go — .9r>, boson + 



1 



2C_1 5D,for 

(7 J) is the D-dimcnsional Stefan-Boltzmann constant^ 



(3.10) 



(JD = 



{D-2){2n)D 



-T{D)aD), 



(3.11) 



where T{D) and C(^) ^-re gamma and zeta functions, 
respectively. Acs^d is the {D — 2)-dimensional effective 
surface area 



(3.12) 



where ftD-2 is the area of the {D — 2)-dimensional 
unit sphere and rg.efi,D is the effective ZJ-dimensional 
Schwarzschild radius in Eq. H3.3(l (with D explicitly 
shown) . 

In the braneworld, one should consider degrees of free- 
dom on the brane and in the bulk separately because 
the former takes 3D phase space and 2D black hole sur- 
face area, while the latter takes 4D phase space and 3D 
surface area. Thus the mass shedding formula takes the 
form below: 



— — ymattcr^cff,4^H ~ 9\ 



graviton 



^cff,5^H 



evap 



for rs < /. (3.13) 



This definition gives half of traditional four-dimensional Stefan- 



Boltzmann constant because 



i?photo 



: 2 was extracted. 



dom of ordinary matter, such as photons, neutrinos, 
electrons, quarks, and their antiparticles etc. Even in 
the braneworld, it should take the same value as 4D 
one. On the other hand, the graviton polarization de- 
gree of freedom reflects the dimensionality as ^graviton = 
[D{D — 3)]/2. We take ^graviton = 5 for 5D black holes. 

What has to be carefully treated is the temperature- 
dependence of the effective degrees of freedom. An in- 
creasing effective degree of freedom of a worming PBH 
will shorten the lifetime of itself. This issue was investi- 
gated by MacGibbon and Webber 0|, and MacGibbon 
|70l |. However, we can expect that the effects on our dis- 
cussion will be limited because our main target is long-life 
PBHs and they have in general lower temperature than 
pion mass around 135 MeV or the so-called quark-hadron 
transition temperature. 

Based on the above arguments, the mass-loss rate for 
a five-dimensional black hole is given as 



/dM 

[~dr 



g: 



2 \lj \ Mi J U 



for M < fMc. (3.14) 
Here the value of the effective degree of freedom gcS,5 is 
J_ fj_ 9C(5) 

2.6 1^160^"'^+ 32^4 5graviton^ 

0.023 for only massless. 



ffeff,5 



The extra factor 1/2.6 in front of the algebraic factors 
is a consequence of numerical calculation carried out by 
Page in the 4D case [l3|- For later use, we also compute 
the 4D effective degrees of freedom, which is defined in 
the 4D evaporation equation as 



dMjt) _ geff.4 (M{t)\ ^ Mi 

dt 3 V Mi J ti ■ 



(3.15) 



4D standard counting for massless species gives goff,4 ~ 
0.00075. 

The solution for Eq. (|3.14() with an appropriate initial 
condition is easily obtained as 



M{t) 



Mp 
Mi 



5cff,5 



t - u 
^7" 



1/2 



Mi 



for M < fM^. (3.16) 



From this formula we can estimate the lifetime of black 
hole tiife. 



^lifc 

4D mass profile is 



5eff,5 \ 



Mi 



ti +tc. 



M{t) 



Mi 



5cff,4 



IT 



1 1/3 



Mi 



(3.17) 



(3.18) 
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As we will see later, all our attention is concentrated 
on the cases in which there are small five-dimensional 
PBHs with sufficiently longer lifetimes than the present 
age of the Universe. For a typical 5D black hole with 
lifetime = to = 13.7 Gyr, its primordial mass M* is 
calculated as 



I 



1031/. 



-1/2 



Thus the corresponding temperature is 



-1/4 



keV. 



(3.19) 



(3.20) 



For the consistency of the above discussion, there is an- 
other requirement that at least a black hole with M* is 
indeed five-dimensional, i.e., rs[M*] < I. It holds when 
I > 1020/4. With / in the typical range 10^0/4-1031/4 of 
our interest, TjJ is 



57 keV < T]5 < 32 MeV, 



(3.21) 



hence the consistency of low temperature approximation 
for those long-life 5D PBHs is confirmed. 



D. The PBH mass spectrum 

In order to evaluate the amount of PBHs in the Uni- 
verse, we take a usual notation iTll 



at (Mi) 



PPBH,Afi jt) 
Prad(t) 



(3.22) 



where ppbh.Mi is the mass density of PBHs with initial 
mass on the order of Mj and prad is the radiation energy 
density. In this paper we only use its primordial value, 
i.e., 

tj,.r\ PPBH,Mi(^i) /n,r\ ^(^0 fo o-l\ 

ai(Mi) = ati{Mi) = — — = at{Mi)——, (3.23) 



Prad(ii) 



ait) 



where ti is given for each A/j from Eq. (|3.1() . The upper 
limit on a; will be written as LIMj. 

Using the above definitions, the comoving number den- 
sity of PBHs with initial mass on the order of Afj is given 

as 

nPBH(M) = ai(M)^^^a(ii)'^^(tiife - to), (3.24) 

where t\i{c is the corresponding lifetime and d{t) is the 
theta function. The relation between npBH(Mi) and the 
differential comoving mass spectrum dn/dMi can be writ- 
ten as mi 



npBaiMi) ~ Mi 



dn 
dMi' 



(3.25) 



Some constraint on npBn{Mi) leads to an equivalent con- 
straint on the mass spectrum dn/dAfj. Thus throughout 
this paper, we use the notion of number density npBH 
and mass spectrum dn/dM in parallel. All the discus- 
sions can be freely converted by the relation 



dn _ TipBH 
dM ^ M ' 



(3.26) 



The relevant mass spectrum for antiprotons, which is 
the issue of this paper, is that at the present moment of 
the Universe. As is mentioned in Appendix IA 31 there are 
four possible destinies for primordial black holes accord- 
ing to their primordial size. In the large extra-dimension 
case / > lO'^o^^^ the present spectrum is 



f-cq 

X /-3/8+^M 



M 



2 , QcbMo 



I 



-25/16- J='/2 



for AC > Af 



^^«i[Afi(A//)]/i/2^/A^ 



5cff,5 



(Af 2 ~ A/2) + Af3 + g,ff,4t0 



-3/2 



5off,5 

for Afc >M> M* 
— ^ai[A/i(Af)]/i/2Af2 [M^+g,e,4to] 
for M > Mr, 



3 < 



-3/2 

(3.27) 



where 



9F 



(tt-F) 



and M: 



I 



(3.28) 



On the other hand in the small extra-dimension case / < 
10^0/4, those lightest PBHs with mass M < M* have 
already died out. 



IV. 



ANTIPROTONS 



In this section, we describe the generation and prop- 
agation of antiprotons in the Galaxy. As our formula- 
tions are solely owed to previous works by some authors 
El El m Hi El m, we here only briefly look back 
the results. Some more details are in Appendix IbI Al- 
most all the differences from the ordinary 4D scenarios 
are Hawking temperature of the PBHs and their mass 
spectrum. 
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Sources 



defined in the following manner: 



1. Antiprotons from PBHs in the galactic halo 

Now let us consider primary antiprotons emitted from 
PBHs distributed over the Milky- Way halo. Hereafter 
we assume spherically symmetric density profile of PBHs 
as dark matter. Based on computational studies on the 
evolution of dark matter density distribution, the present 
profile is parameterized in some ways. One of its general 
expressions in the cylindrical coordinates (r = is the 
galactic center and z = corresponds to the disc) is 



pPBH(i?0,O) 
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(4.1) 
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where subscript j = q,q, g, stands for quark or 

gauge boson species emitted from a black hole. dg^x/dE 
is the differential fragmentation ratio into antihadrons 
for a particle j. We included antineutrons because they 
immediately decay into antiprotons. 

From above considerations, we are now able to write 
down the differential antiproton emission rate from pri- 
mordial black holes as 



gP"(r,z,i?) 



'^,(M,E)§ir,z,M)dM. (4.4) 



where i?corc is the so-called core radius and Rq is the dis- 
tance from the galactic center to the solar system [73 . [73j . 
We here use a = 2, 7 = 0, e = 1, and i?core = 3 kpc. 
Choice of parameters is expected not to significantly af- 
fect the results much '48^. Which mass spectrum should 
be used is determined by how long it takes for the mes- 
sages to propagate from evaporating PBHs up to the in- 
strument on the earth. For antiprotons, the time scale 
is sufficiently shorter than the cosmic age, hence we take 
the present mass spectrum dn/dM given in Eq. (|3.27() . A 
mass spectrum has such a value that is averaged over the 
whole Universe, so now let us define the effective mass 
spectrum in the Galaxy dn/dAf as 



dh 
dM 
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(4.2) 



Typical primary emissivity at the solar neighborhood is 
shown in Fig. ^ 




1.0 

Kinetic energy [GeV] 



where G represents the enhancement of the galactic mat- 
ter density measured at the solar neighborhood relative 
to the extragalactic value. It is said that Q ^ 10^ ^M- 

Next, we consider particle emission degrees of freedom 
from black holes of quarks, antiquarks, leptons, antilep- 
tons, and gauge bosons. Among them, quarks, anti- 
quarks, and gauge bosons are able to generate antipro- 
tons via hadronization process in quark or gluon jets. We 
make use of pythia/jetset jT^ for it. Our procedure is 
not so different from the 4D situation; the necessary in- 
gredients that we must fix are only the Hawking tempera- 
ture and Schwarzschild radius of a black hole. Their only 
difference from the four-dimensional case is 0(1) numeri- 
cal factor coming from the topological nature of the black 
hole space-time. Hence in the calculation, we can trust- 
fully use the 4D result of spin and energy-dependence of 
the cross section, which was explicitly shown by MacGib- 
bon and Webber [? ] . From here we denote the effective 
differential emission rate of antiprotons as d^Np/dEdt 



FIG. 1: Illustration of primary emissivity of antiprotons. 
The difference between the two is mass spectrum of PBH. 
Dotted line is 5D one in a typical case I = lO^^h, F = 1.0 
and Gai = 10"^^ 



2. Secondary antiprotons as the background 

Next we consider secondary antiprotons. The sec- 
ondary antiprotons are generated in the particle reaction 
within the galactic disk between cosmic rays and the nu- 
clei in the interstellar gas, such as N -I- N' ^ p + X . 
These process is thought to be well described by the 
known experimental data, so we treat secondaries as the 
background. Actually this secondary component domi- 
nates the observed antiproton flux. Because of the cos- 
mic abundances, the most dominant sources are collisions 



between cosmic-ray protons and 90% hydrogen plus 10% 
helium in the interstellar gas [s^Ff^ : 

pR ^ pX, 
pHe — > pX. 

Proton cosmic ray has been measured by observations 
such as IMAX 0LBESS [zi, and Alpha Magnetic Spec- 
trometer (AMS) [23 ■ The result of the observation of 
high-energy cosmic-ray protons at the solar neighbor- 
hood is well fitted by [S^l 



g*"(r,0,i5p) 
= (riH + 42/3nHc) 



_non — anil /rp/ 



[4^<I>p(r,0,4)] d^ 



- (nn + 42/3nHc)a^7-'^""(Tp) [4^$p(r, 0, E^)] . (4.8) 

Due to the weakness of those interactions, the tertiary 
source is left ignored until numerical analyses. 



N { ^J" \ m~2 sr-i s"i GeV" 



1 GeV 



(4.5) 

with parameterizations N — 13249 and 7 = 2.72, where 
Tp = Ep — Trip is kinetic energy of protons. Here we 
assume it is constant for r and use the above value. 

As for the cross section, we make use of the most recent 
parameterization made by Duperray et al. [80l |. which 
gives better than the well-known result by Tan and 
The differential cross section is given as 



dUpA- 



dEp 



= 2 / 2'KppSiii9p [ Ep—^ 



dOp, 



(4.6) 



where pp — y Ep — nip and 9p are momentum and az- 

imuthal angle in the laboratory frame, respectively. The 
limits of 9p are determined by kinematics. Again we mul- 
tiplied the expression by 2 to incorporate antineutrons. 

Finally, we obtain the production rate of the secondary 
antiprotons in the galactic disc as 



00 

q^'^^r, 0, Ep) = J dEp [47r$p(r, Ep)] 
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X ^ HA 

A=H,Ho 



dapA^px 
dEf, 



iEp,Ep), (4.7) 



where the threshold of the interactions Eth is Snip. We 
use riH = 1 cm^'^ and mjn = 0.1 cm~^ for number the 
densities of the targets |7J| . 



3. The tertiary source 

There is another source of antiprotons, so-called ter- 
tiary component. This is not a real production mecha- 
nism of antiproton but just energy exchanges between 
propagating antiprotons and ambient protons in the 
galactic disc. The relevant interaction is the resonant 
excitation of the proton, which is non-annihilating but 
inelastic ^S^. The tertiary component is expressed as 



B. Analyses on antiproton flux 

It is obvious that p's are only emitted from PBHs with 
typical temperature Th ^ rrip ^ 1 GeV. The mass corre- 
sponding with Th 1 GeV, MgcV, is from Eq. ^^THl 



JGoV 
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g for / > 10^" (4.9) 
for / < 10^". (4.10) 



A condition MqoV < -^^p holds in any situation, which 
means that PBHs being such hot by nature, i.e., since 
their formation, cannot survive till now. Thus contri- 
butions to the antiproton flux are all from evaporating 
region in the present mass spectrum, and we further find 
that the contribution to antiproton flux is dominated by 
only those PBHs with mass around MceV- As a conse- 
quence, we understand that no braneworld signature in 
the antiproton flux is expected when I < 10^" because 
then the shape of the present mass spectrum around 
MceV does not deviate at all from that in the 4D case. 
On the other hand, when I > 10^°, MqcV PBHs are 
five dimensional and the effect of the braneworld will be 
strongly reflected in the p flux. 

Combining with another fact that the p emission rate 
per a PBH is completely determined by its Hawking tem- 
perature only, one reaches the conclusion: the spectral 
shape, such as peak location of the resultant antiproton 
flux, is never changed even if ^ or varies, while the 
only possible modification is over-all fluctuation. Under 
the low-mass limit MceV M*, we can schematically 
write down the dependence of the observable p flux at 
the solar neighborhood, ^p, on the mass spectrum as 
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d^N dfi 
dEdtdM 

d'^N M 
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(4.11) 
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with the mass spectrum from Eq. (|A18|I : 
3 



mum value Zmin- Then the above proportionahty becomes 



dn 
dM 



Mgov 
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(4.12) 



where MqoV c< ^ ^- Finally, the parameter dependence 
of the flux is obtained as 



(4.13) 



Here we show the evolution of PBH mass spectrum in a 
typical case I = 10^^Z4 and F — 1.0 as in Fig. [5] 



5D primordial 
5D present 
4D present 




10" 10" 10' 
PBH mass [g] 



FIG. 2: Lines are 5D PBH mass spectrum soon after accre- 
tion (solid), 5D present one (dotted), and present one in the 
4D case (dashed), respectively. Braneworld parameters are 
taken to be / = 10^^/4 and F = 1.0. 



The most important result is that in general the an- 
tiproton flux is a decreasing function of I. Recalling 
that F ~ %F/{^Ti — 8F) takes a value in the range 
of < < 0.53, it is understood that the index of 
I, which is denoted as p, is in fact always negative: 
p e [-0.8125,-0.0175]. Hence we find that abundant 
PBHs set a lower bound on the bulk curvature radius 
I. In another word, the limits on ai is proportional to 
^i3/i6-3.?'/2^ which is an increasing function of I. Large I 
can relax the upper bound on ai. 

For completeness, we examine how the F dependence 
is. Taking a logarithm of Eq. (|4.13|) . we obtain 

/13 3F\ F 
log $p oc log [Gai] - ( Ye " ~ j ^ ~ "2 [gcsMo] ■ 

(4.14) 

When the flux reaches some upper bound, which means 
that cki also reaches its upper limit LIMi, I has the mini- 



log [U.L.] = log [gUM: 
13 3F 
16 ~ ~ 



log /„ 



F 



log [goff,5io] , (4.15) 



where U.L. is an undetermined constant which corre- 
sponds with the upper limit flux. In principle U.L. cannot 
be determined unless the propagation process is solved. 
However, even without the knowledge of propagation, we 
can find more about the above relation. What is impor- 
tant is the fact that the only difference for each case is the 
over-all fluctuation of the antiproton flux. Therefore we 
can understand that U.L. is a universal value. Keeping 
it in mind, we modify the relation as 



log^i, 



_ log[.9cff,5^o] 
3 

481og[gLIMi/U.L.]/log[geff.5^o] - 13 
13 - 24J^ 



(4.16) 



Although wc cannot determine in general the sign of the 
numerator in the right hand side, we are particularly 
interested in the cases of large PBH abundance, which 
means that the numerator is positive. For the same rea- 
son as the flux is a decreasing function of the denomi- 
nator is always positive for < F < 0.53. Therefore we 
finally find that the lower bounds for I draw a family of 
hyperbolas expressed in Eq. H4.16|) . 



C. Numerical studies 

Our final task is to give quantified conclusions by a nu- 
merical study on antiproton propagation in the Galaxy. 
We adopt here a semi-analytical method which was first 
proposed by Webber, Lee, and Gupta 1501 and has been 
extensively utilized by other authors |48l l5lL Is^ . IsM Issl 



1. Diffusion 

This method is based on the diffusion model which was 
proposed by Ginzburg, Khazan, and Ptuskin [49j . 

The full diffusion equation in the steady-state for the 
differential number density N(r,z,E) = dup/dE with 
the source Q is js^ 

dN ^ ^ 
= — = V • [K{E)VN{r, z, E) - V,{r, z)N{r, z, E)] 
at 



v-K(r,z) d 



dE 
d 
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Nir,z,E) 



T{E)N{r,z,E) 



Q{r, z,E) + — {- [b,,,cciE) + bcon{E)] N{r, z, E)} 



d_ 
'dE 



f}'Kpp{E)^N{r,z,E) 



(4.17) 
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K and Kpp are the diffusion coefficient in the real space 
and energy space, respectively. We assume power-law for 
the coefficient as K = KqPp^ with 5 = 0.6 |E1 El El 
El El 113. In terms of K, we put Kpp (x V^p'^ /K, where 
Va is the Alfven velocity. Vc is the velocity of convective 
wind, r expresses the escape probability from an energy 
bin via general inelastic interactions, while 6rcadd and 
^coii correspond with energy losses by reacceleration and 
collision, respectively. 

We put a cylindrically symmetric region, where dif- 
fusion occurs. It is separated into two parts: the thin 
galactic disc, and the so-called diffusion halo sandwich- 
ing the disc. The radius of the cylinder, R, is set to 
20 kpc and the solar system is located at Rq = 8 kpc 
from the galactic center [s^ . Issf . The half- height of the 
disc, denoted as h, is fixed at 100 pc, while that of the 
diffusion halo, L, is set to 6 kpc. Owing to the cylindri- 
cal symmetry, we express all the physical quantities that 
have radial dependence in terms of the Bessel functions 
of order zero, so that the diffusion equation is reduced to 
those for the Bessel components. After appropriate eval- 
uation of the terms, which is presented in Appendix IbI 
the final form of the diffusion equation for the antiproton 
number density becomes 



= 



N,{z,E) 



- V 

K{E)^ + 2h5{z)Tpp{E) 



N,{z,E) 



-biE)+p'Kpp{E)^ 



N,{z,E). 



(z = l,2,...) (4.18) 

b = 6reacc + &coii + &adiab, whcrc 5adiab correspouds with 
energy loss by adiabatic expansion. This is the basic 
equation to be solved. 

2. Solar modulation 

An obstacle for antiprotons approaching the earth is 
the solar wind. The solar wind is a mixture of plasma 
and frozen-in magnetic field spouting from the sun. As a 
result of propagation against the wind, a significant part 
of antiproton energy will be lost and the flux measured 
at the top of atmosphere (TOA) on the earth will be 
suppressed relative to the value at interstellar (IS) space. 
This phenomenon is called solar modulation, which is 
synchronized with the solar activity with 11 years peri- 
odicity. The most energetic period of the sun is called 
solar maximum, while the least energetic is solar mini- 
mum. 

Although there are several ways to treat solar modula- 
tion, it has been known that they have their merits and 
demerits |8^, so we here use a basic method, so-called 
force field approximation. Basically, in every framework 
for the process, an electrical-potential-like quantity 0, 



which has a unit of GV. is intoduced to express the dif- 
ficulty for cosmic-rays to flow against the solar wind. 
Actually the force-field approximation method does. In 
practice, (j) is estimated by the observations of cosmic- 
ray proton or neutron flux. This parameter takes a value 
from 500 MV at the solar minimum to 1 GV at the solar 
maximum. In the model of force-field approximation, the 
resultant energy at TOA is determined by that at IS and 
6 as 



^TOA 

A 



E^ 
A 



\Ze\4> 
A ' 



(4.19) 



where A and Z are mass and the charge number of the 
cosmic-ray content; for antiprotons they are both set to 
unity. The ratio of the flux is written in terms of their 
momentum as 



^TOA /pTOA 



$is 



p 



IS 



(4.20) 



The solar modulation is so crucial for our discussion 
because it generally results in the energy loss ranging 
from a few hundred MeV at solar minimum to over GeV 
at solar maximum. Hence most signatures of PBH evap- 
oration, typically appearing in sub-GeV region, will be 
swept out at solar maximum. Therefore we should take 
solar minimum data to be compared with the theoreti- 
cal framework. This is the reason why we will concen- 
trate on BESS and CAPRICE data at the solar minimum 

mEllillillii. 



3. Semi- analytical method 

The diffusion equation can be regarded as a composite 
of the two qualitatively distinctive parts, i.e., diffusion in 
the real space and that in the energy space. Basically, 
the energy part can be treated as a perturbation. 

To integrate such a type of differential equation, we 
use a semi-analytical method [s^- First we integrate 
the equation in the direction of z, omitting both energy 
gain/loss terms and the tertiary source. We use the an- 
nihilation cross section for Vpp here (see Appendix IB 2} . 
Then that solution is substituted into resulting pertur- 
bative differential equation in the energy space together 
with the tertiary source term, and we solve it numeri- 
cally^. We will iterate the energy diffusion till the solu- 
tion converges. 



4- Choice of propagation parameters 

In order to clarify the issue, we shall fix values for the 
free parameters: three diffusion parameters and one so- 



^ In this final process, we substitute the total cross section for 
inelastic interaction into Fao. 
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lar modulation parameter, i.e., Kq^ Va, Vc, and (j). We 
understand that the parameters should be in principle 
determined independently of the braneworld. Moreover, 
our purpose is solely focused on seeing how brane pa- 
rameters affect the antiproton flux. It means that mod- 
els for propagation is here only required to be reasonable, 
thus we do not have to repeat fluctuating propagation pa- 
rameters and searching the least for every braneworld 
parameter set. It is sufficient for our purpose to fix rea- 
sonable values for the above parameters. 

First we consider diffusion. The three quantities should 
be in principle determined by observations, but unfor- 
tunately there has been no direct, precise observational 
result. Thus the only trustful wa y to fix them is fitting 
cosmic-ray data to the theory |4M|52j. In the procedure, 
basically we must try to exclude uncertainties originat- 
ing from cosmic-ray sources otherwise they will introduce 
new errors for the parameters. From such a point of view, 
antiproton is not a very appropriate cosmic-ray nucleus 
because we have been taking a standpoint such that their 
production rate is largely varied by both the amount of 
PBHs and braneworld parameters. Therefore we borrow 
a result obtained by the analyses on the cosmic-ray car- 
bon flux to boron (B/C) ratio [Hilil. Although the 
authors could not reach the goal to get the best-fit pa- 
rameter set, they succeeded to show sufficiently reason- 
able parameter regions. Our choice for the parameters, 
which is shown in Tab.Ql is consistent with their result. 

The next problem is the solar modulation. As was 
mentioned before, we use only solar minimum data within 
a few years. Although each author having performed 
every year's observation has determined the best-fit 4> 
for the solar activity, their method to model the solar 
modulation is respectively different; especially, some are 
not the same as our choice, the force-field approximation. 
Hence we cannot directly use their values. Here we take 
a conservative value 500 MV for (h. 



TABLE I: Propagation parameters sufficiently good for our 
purpose. They are all consistent with previous results ob- 



tained for ordinary cosmic rays |53l l55ll5i 



Quantity 


Ko Va 


Vc 


4> 


Value 


0.012 kpc^ Myr"^ 60 km s"^ 


15 km s" 


500 MV 



BESS 1995 I — H 

BESS 1997 ^-y. 

BESS 1998 i 

CAPRICE1994 : e 

CAPRICE 1998 t- -m 

total 

secondary 

-I primary 

primary (4D) 




1.0 10.0 
Kinetic energy [GeV] 



FIG. 3: Total antiproton flux with a typical braneworld 
parameter set I = 10'^^ Z4 and F — 1.0. Primary flux is 
the best-fit one {gai = 1.1 X 10~^^). 4D best-fit primary 
(Qai = 1.8 X 10"^'*) is also shown for comparison (4D total is 
not displayed). 



From Fig. 01 we find that a few percent precision ob- 
servation in sub-GeV region can distinguish the two pri- 
mary flux, i.e., we will be able to test the braneworld 
with future observations. 

The final results of this section are allowed regions in 
the parameter plane {l,F) for best-fit (Fig.g)), 50% C.L. 
(Fig.EJ), and 90% C.L. (Fig.©. 
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5. Allowed regions 

The final project is to obtain upper limits on the abun- 
dance of PBHs for each braneworld parameter set (l,F). 
We fix the diffusion parameters and only let the PBH 
abundance Qa-i move. First we show spectra with brane 
parameters I = IQ^^l^ and F = 1.0 in Fig. |21 which are 
corresponding to the best-fit Qai. 



FIG. 4: Boundaries of [l, F) in the braneworld-parameter 
plane for PBH abundance go.; = 10"^^-10"^^ (from top to 
bottom) using best-fit parameter. They are all lower bounds 
of I, i.e., left-top regions are allowed for each PBH abundance. 

At this stage, we can make a crude estimation of 
the cosmological PBH abundance. We see in Figs. 0]- 
Elthat non-trivial results are obtained for galactic PBH 
abundance ranging in the region 10~^^ > Qa^ > 10"^^. 
Although Q is not determined yet, its order of magni- 
tude is around 10^, so the results are indicating that 
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FIG. 5; Boundaries of {l,F) for Qai = 10~"-10"^^ 
top to bottom) using 50% C.L. data. 
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FIG. 



Boundaries of {l,F) for Gai = 10""-10"^^ (from 



top to bottom) using 90% C.L. data. 



LIMi - 10^1^-10-28 for I > 1020^4. Nothing happens 
in the case ai < lO^^^. Comparing the constraint with 
that in the diffuse photon case jH3 1 we find that they are 
almost the same order constraints. 

We here comment on the result obtained in the 4D 
setup. The best-fit PBH abundance is Qa-i « 1.8 x 10^^"', 
which is consistent with the result obtained by Barrau et 
al. When we assume Q ^ 10^, ai becomes around 

10"^^. Although there are large uncertainties, it is about 
two orders of magnitude lower than the upper bound ob- 
tained from the diffuse gamma-ray background |^ . 



V. DISCUSSIONS AND CONCLUSIONS 

After the analyses, free parameters /, F, and ai are 
left undetermined. They both originate from the brane 
early Universe but can be in principle independently de- 



termined; the former is a problem of gravitational col- 
lapse and the latter is that of relativistic hydrodynam- 
ics on the gravitational background, ctj is directly con- 
nected with the primordial density perturbation spectra, 
and can be potentially given under specific braneworld 
models, such as the bulk inflaton model [s^l or the cyclic 
Universe model and so on. Note that our results will 
not be qualitatively changed even if the assumption of 
constant ai is moderately violated; in the both analyses, 
what was important was only an absolute value of the 
mass spectrum for narrow mass range. 

In the analyses, we assumed that an evaporating black 
hole holds the same temperature as the initial value 
through its lifetime. Hence the effective total degree of 
freedom, ^eff, was fixed at a constant, which is for zero 
temperature. However, as was shown by MacGibbon [70l |. 
the degree increases maximally by an order of magnitude. 
Due to such an underestimate of evaporation rate of hot 
PBHs, our resultant antiproton flux was somewhat over- 
estimated and constraints were in turn virtually strength- 
ened. 

We should make a remark on the range of viability of 
our methodology presented in this and the previous pa- 
per [355 which served constraints on a higher-dimensional 
cosmology. Although our interest has been throughout 
concentrated on a particular cosmological model of RS2 
braneworld, we insist that the use of Hawking radiation 
from small black holes is available for a broad range 
of models with large extra dimensions. The most cru- 
cial feature for the purpose is the change of Hawking- 
radiation spectra emerging when a hole notices extra 
space extending around. This will generally occurs if 
only the target theory has compact extra dimensions 
and the hole shrinks enough below the compactifica- 
tion scale. This compactification is not necessarily re- 
alized by periodic boundary conditions, but RS2-like ef- 
fective localization of gravity is also permissible. Exam- 
ples of phenomenological models contain Arkani-Hamed- 
Dimopoulos-Dvali 89] and Hoi^ava-Witten cosmologies 
[90| . Further, if the theory has a fundamental mass scale 
lower than 4D Planck, which is somewhat likely, then 
Hawking temperature will drops like that in the RS2 
analysis and hence our results may be directly applied 
with some modification. 

Here we summarize this paper. We investigated an- 
tiproton flux emitted from PBHs distributed over the 
Galaxy as dark matter. Mass of PBHs which could con- 
tribute to the antiproton spectrum was in the 5D evapo- 
rating region of the present mass spectrum, which meant 
that the effect of accretion had been almost completely 
swept out. Because the antiproton emissivity of individ- 
ual PBH was almost unchanged from that in the 4D case, 
the change in the resultant antiproton flux only came 
from the change of the slope of the mass spectrum; in 
the large extra-dimension case, dn/dM ~ M, while with 
the small extra-dimension, dn/dM ~ M^. However, this 
also does not affect the antiproton flux so much because 
the flux is still dominated by the contribution of PBHs 
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with temperatures around 1 GeV. Finally we understood 
that the whole situation is almost identical with the 4D 
case, which was shown by the expected flux in Fig. 13 
Therefore we expect that even future high-precision ex- 
periments are done in sub-GeV kinetic-energy region, it 
will be difficult to identify braneworld signatures from 
the spectral shape of antiproton flux only; this is the 
crucial difference from the photon case obtained in our 
previous paper [35l| . 

As the most impressive result, we found that the an- 
tiproton flux from brane PBHs was a decreasing function 
of the bulk curvature radius I like 



(X «./-13/16+3^/2 



(5.1) 



The reason why the flux decreases is that the temperature 
of a black hole is lowered as I enlarges, while accretion 
takes a quite limited role in contrast to the photon case 
[s^. Therefore lower bounds on the size of the extra 
dimension were obtained. Of course, on the contrary, 
constraints on ai from antiprotons could be relaxed by 
the factor of /i3/i6-3J='/2 ^f^gj. allowed regions of the 
braneworld parameters from antiproton observation were 
presented in Figs. QHH for each confidence level. 
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FIG. 7: The composite of the two results of allowed regions 
obtained in Sec.|l V|and in our previous paper for high-energy 
diffuse photons 35] (presentation is a bit refined). Regions 
allowed by both observations is in the left-hand side of each 
boundary. The density enhancement in the Galaxy Q is set 
to 10^ here. Numbers in the figure indicate Ofi's. For ai > 
10~^®, diffuse photon constraints are stronger than that of 
antiprotons. However, for a\ ~ 10~^®-10~^^, antiprotons set 
lower bounds on the size of the extra dimension, which was 
not done by diffuse photons. For a\ < 10"'^*, antiprotons do 
almost nothing. 

Finally we show the composite of the allowed regions 
obtained from two distinctive observations, i.e., above- 
mentioned antiprotons and the diffuse photon back- 
ground [S^; see Fig. [7| In Sec. IIVI the density en- 
hancement of the Galaxy Q was left undetermined, but 
we here use a conservative value Q ^ 10^. Thus the 



boundaries for antiprotons in Fig. |7| now correspond to 
a; = 10-18-10-28. We find two facts from the fig- 
ure: (i) for boundaries corresponding to ci\ > 10^^^, 
constraints from diffuse photons are always more strict 
than that from antiprotons, while (ii) for boundaries of 
Q!i ^ 10-2^-10-28^ they have chance to cooperate and set 
the lower bound of I. 

In the given situation, the strongest implication is ob- 
tained when we assume that the apparent low-energy en- 
hancement of the antiproton flux is actually due to the 
PBH evaporation. If so, it means that the PBH abun- 
dance tti is near the upper bound, at most one order 
of magnitude below it. On the other hand, we here 
should remind that there is another requirement that 
diffuse photons from PBHs must be hidden in the back- 
ground component. Referring to Fig.[7| such a condition 
is interpreted as follows. Sub-critical PBH abundance 
(tti ~ LIMj) means that allowed braneworld parameters 
are nearly restricted in a belt-shaped region expanding 
slightly above antiproton- lower bounds, i.e., dotted lines. 
On the other hand, the parameters must be in a corre- 
sponding allowed region determined by photons, which is 
on the left-hand side of dashed lines. We find that such 
two conditions are only satisfied when ai « lO-^^-lQ-^^. 
Furthermore, it means I < lO^"*. This upper bound is 
about seven orders of magnitude stronger than that ob- 
tained in today's experiments, and may become a dis- 
advantage of the braneworld scenarios. In the future, 
more contamination-free observations of PBHs, such as 
antideuterons or other anti nuclei, may give an answer to 
the problem. 
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APPENDIX A: SHAPES OF THE PBH MASS 
SPECTRUM 

1. Initial spectrum 

To begin with, we calculate the initial mass spectrum. 
The first case is with mass Mj < fMc- PBHs with such 
light mass are formed in the brane high-energy phase 
and experience the accretion. From Eqs. (|2.17|) . (|3.1|l . 
and H3.24|l . the comoving mass spectrum of PBHs with 
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initial mass Mi is obtained as 
— ai(Mi)- 



dMi 



Mi Mi 

- .^,^.wl/8;-3/8^j.-17/8 

- OlV^TT +3/2 i 

for Mi < /Mc. (Al) 

This has a slightly red spectral index in comparison with 
the 4D case. For fom'-dimensional PBHs with Mi > /Mc, 
from Eqs. (jTT7|) . JTIJ, and 



for Mi > /Mc 



(A2) 



This mass dependence is actually the same as that of the 
pure 4D case. 



3. Present spectrum 

The present mass spectrum is a bit complicated be- 
cause evaporation is taken into account. After accretion 
on PBHs ends, PBHs begin to radiate and lose their mass. 
An important feature is that the mass decreasing rate of 
each PBH is different according to its dimensionality. 

In order to obtain the complete mass spectrum, we ar- 
gue detailed evolution of each PBH. We can understand 
that at the present age of the Universe, there are four 
destinies prepared for each PBH according to its primor- 
dial mass Mpi (i) born as a 4D PBH and now still be 
4D one (their mass is symbolically denoted as Mp), (ii) 
born as a 4D PBH but now is 5D one (Mp), (iii) born as 
a 5D PBH and now evaporating (Mp"), and (iv) born as 
a 5D PBH and having already evaporated away (Mp^). 
The mass hierarchy among above typical values is 



2. Primordial spectrum 

The next task we have to do is to consider the effects of 
accretion [s^, IssL l35l| . Here we call the quantities which 
take accretion into consideration but does not evapora- 
tion as "primordial" ones. 

5D black holes are born before tc and their masses will 
increase due to accretion during the brane high-energy 
era. Resulting primordial mass Mp after accretion is ob- 
tained from Eqs. (|2.17() and H3.1|l as 

2F/TT / r, \ F/lT 



M„ - I ^ 



Mi = 4Af 



AL. 



Combining it with Eq. (|Aip . we obtain the primordial 
PBH mass spectrum as 

dn dn dM 



dMp dMi dMp 



^^«i[A/i(Mp)] (l + ^) 



for Mp<fM„ 



(A4) 

9i^/(87r - 8F) 



where a new accretion parameter T 
was introduced for convenience. Particularly we find a 
relation between the distorted spectrum and that with 
no accretion as follows: 

2 /; \ ^ dn 



dn 



a; (Mi) 



dMn 



4A// 



M, 



dMn 



(A5) 



no accretion 



Ct;i(A/p) y jvjp y ^^jvjp 

In the A/j > /Mc case, they are born in the standard 
radiation-dominated era and do not experience accretion, 
so this part of the primordial mass spectrum is unchanged 
from the initial spectrum. Hence, identifying Afp = Afj, 

^ = ^ for A/p>/Mc, (A6) 
dMp dAfi P - ./ c, K J 

where dn/dM{ is given in Eq. ijUJ. 



A/;, > AP^ > Af- 



> M'' 



(A7) 



What is non-trivial is the relations between every typical 
mass and M* , the primordial mass of a PBH with lifetime 
equal to the age of the Universe to. If a condition M* > 
Mp' holds, PBHs which were originally five-dimensional 
have already died out and there are only originally four- 
dimensional ones now. To the contrary, in the case of 
Af* < M'^\ there are still originally 5D PBHs, which 
have experienced accretion in the brane early Universe. 

Here we note that our framework, braneworld, itself 
has a mass scale Mc, which determines whether a BH is 
five dimensional or four dimensional according to M ^ 
Afc. If we require that PBHs having just evaporated were 
originally five-dimensional, we need Mc > M* and hence 
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(A8) 



This condition is in general important for constraints 
from cosmic-rays since, if it does not hold, we cannot 
see any braneworld signature in those observations; see 
Sec. lIVI fthere a bit stringent condition was required) and 
our previous paper 35]. Additionally, there is another 
mass scale implemented in the framework. A/*, which is 
defined as 



m: = A M2 



I ■ 



(A9) 



Recalling Eq. (|3.1t)|) . it is understood that a PBH 
presently with M* has primordial mass Mp = Mc. If 
the mass of a PBH is now M < M* , we can see that it 
was born as 5D one because the condition means that its 
primordial mass Afp was smaller than M^. 

Now let us investigate the resultant mass spectrum. 
We already understood that PBHs in the mass range 
smaller than M* are all on the way to evaporation. Com- 
bining the expression of the primordial mass spectrum 
Eq. (jA4|) and the mass evolution formula Eq. H3.16|l . we 



15 



can calculate the present mass spectrum, which is for- 
mally written as 



dn 
dM 



dn dMi dMp 
dMi dAfp dM 



(AlO) 



For the three types of PBH (i)-(iii), we can calculate 
derivative factors in the above equation and finally obtain 
whole mass spectrum as follows. Note that we fix / = 1 
for simplicity. 



a. M>M^ 

The present mass of the black hole is larger than 
boundary value. It is now four-dimensional and was born 
as 4D one. Its history is described as 

Mp = Mi and M{t) = [Af^ - g,s,^{t - t,)] ^'^ . 

(All) 

Thus the present mass spectrum of this region becomes 

to" 



i-cq 

for M > Mr 



(A12) 



6. Mr > M > M* 



Next is the most complicated case; large initial mass, 
but not enough. A PBH in this category was born as a 
4D one and now five-dimcnsional. Such ones only exist 
if I ~ 10^^, so usually we do not have to pay attention to 
them. The history is 



and 



M{t) 



ffoff,5 



Mp = Mi 



{Ml - M3) - t 



M^ 



The present mass spectrum is 



(A13) 



1/2 



(A14) 



dM 167r fij^ 5cff,5 



1 On 



X M 



9cSA 



I (M^ - M2) + Af3 + grs^4, 



for Mc > M > M* 



-3/2 



(A15) 



of PBH can only exist under a particular condition on the 
bulk radius, Eq. (|A8|I . The history is 



Mp = {4:fiy'''Ml 



and 



M{t) = 



p I 

Hence the present mass spectrum is 



1/2 



(A16) 
(A17) 



dn 



''CQ 



ai[Mi(Afp)] 1 



for M* > M. 



M2 + 



-25/16-:p/2 

(A18) 



APPENDIX B: PROPAGATION OF 
ANTIPROTONS IN THE GALAXY 

In this paper, we adopt a semi-analytical method which 
was first proposed by Webber, Lee, and Gupta 's^ and 
has been extensively utilized by other authors ;48t . .51i..53 . 
Is^ Issl Is^ . This method is based on the diffusion model 
proposed by Ginzburg, Khazan, and Ptuskin '49'|. 

At the beginning, we write down the full diffusion equa- 
tion in the steady-state for the differential number den- 
sity N{r, z, E) = dnp/dE 

dN ^ ^ 
= — = V • [K{E)VN{r, z, E) ~ V,{r, z)N{r, z, E)] 



dE 

d 



^N{r,z,E) 



T{E)N{r,z,E) 



+ Q{r,z,E) + — {-[b 

rcacc 

{E)+broix{E)]N{r,z,E)} 



d_ 

'dE 



P^KppiE)-^N{r,z,E) 



(Bl) 



As was mentioned in Sec. IIVI we assume cylindrical 
symmetry for the Galaxy, and expand all the physical 
quantities in terms of the Bessel functions of order zero. 
Such quantities, for example, differential number density 
N{r, z, E) or the source term Q(r, z, E), can be expanded 
as 

oo 

/(r,z,i?) = ^/,(z,i?)Jo(c.^), (B2) 

i—l 

where the components are given by 
R 

2 Ft / t \ ciT 



c. M* > M 

A PBH of this type is sufficiently light and was natively 
five-dimensional. Readers should remind that such a kind 



and (i 's are the zeros of Jq. As a consequence, the prob- 
lem has been reduced to diffusion along the vertical di- 
rection z. From now, let us see the roles of each term in 
the diffusion equation Eq. (|B1|) . 



16 



Diffusion and convection 



3. Energy gain — reacceleration 



K{E) and Vc^ which appeared in the first and sec- 
ond term in the diffusion equation, stand for diffusion 
and convection, respectively. During the propagation in 
the halo, charged particles are diffused by the random 
galactic magnetic field. We adopt standard discussions 
[Hi m m iHl m m 113 for it; the energy-dependent 
diffusion coefficient is given as 



K{E) = Kopn^ 



(B4) 



where /3 denotes the velocity and TZ = p/ Z is the rigidity. 
For antiprotons, TZ is identical with its momentum p. 
Here, in principle is only an undetermined parameter 
since 5 should be fixed by microphysics. We have set 5 
to 0.6 throughout this pap er in accordance with previous 

results m m m H iHi m . 

Another potentially important effect for cosmic-ray 
transfer is convection. The convective (or galactic) wind, 
which consists of plasma medium, flows outward from the 
galactic disc 92j . For simplicity, we assume that it holds 
a constant value in the whole region of the diffusion halo, 
i.e., the convective wind Vc takes the form 



(B5) 



where T4 is a parameter to be fixed. The convective 
flow also contributes to the energy loss in the third term 
through the adiabatic expansion of the plasma. This en- 
ergy loss takes place only in the region where its diver- 
gence has a non-zero value, namely on the galactic disc. 
The adiabatic energy-loss expressed by the third term is 
eventually written down as 

6adiab = -2M(z)T(^^±^)^, (B6) 

where T and m are kinetic energy and rest mass of an- 
tiprotons, respectively. 



2. Inelastic interactions 

The fourth term represents the total rate of interac- 
tions with the interstellar gas by which an antiproton es- 
capes from the energy bin. Such interactions occur only 
in the galactic disc. For collisions between antiprotons 
and hydrogen nuclei, r(i?) is defined as 

T{E) =. 2hS{z)Tpp{E), (B7) 



where 



Tpp{E) = (nn + A''^^nuo)p(rpp{E). (B8) 



app is the total cross section of the interaction of interest. 



During propagation, antiprotons are reaccelerated by 
the random magnetic field. The effect of the Fermi ac- 
celeration appears throu gh t he diffusion coefficient in the 
energy space defined as jSJ] 



Kpp{E) = 



3SiA-5^){A-S) K{E)'' 



(B9) 



where Va is the Alfven velocity of the galactic magnetic 
field. 

Two terms in the diffusion equation contain Kpp: one 
is the coefficient of the first derivative of energy 



1 + 13^ 

^rcacc = 2hS{z) — Kpp{E), 

and another is in the second derivative of energy. 



(BIO) 



4. Energy losses 

We take into account three types of energy- loss for 6coii: 
two collisional processes, namely Coulomb scattering and 
ionization, and one adiabatic process. For these energy 
losses, see 93, 94J. 



5. Sources 

Three kinds of the sources have been described in 
Sec. lIVI With the thin-disc approximation on the Galaxy, 
the source term takes a form as 

= qf + 2/i<5(z)qr + 2h5{z)qrm, (BU) 



where 

?r'(=.E) 



= (0,i?)= / dEpPNp,,{Ep) J2 

A U 1 



A=H,He 



dE 



(B12) 

{Ep , E) , 
(B13) 



and 

qrmio,E) 

.2/3 



non — aim ('j^f\ 

nH + 4^/^nHc) / dE' ^, ^ ' m{z,E') 



T' 



(nn -f42/3nHc)CT" 



non— aim 
pp 



{T)PN,{z,E). (B14) 



dfii/dM's in Eq. (|B12p are Bessel components oidn/dM. 
A relation between flux and number density 



<i>{r,z,E) = ^N{r,z,E) 



(B15) 



was made use of. 
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6. Final form of the diffusion equation 

After above considerations, the final form of the diffu- 
sion equation for the antiproton number density becomes 



where b = fercacc + ''coii + &adiab- This is the basic equation 
to be solved. 



d 



Ni{z,E) 



K{E)^ + 2h5{z)Vpj,{E) 



-h{E) + P''Kpp{E) 



Ni{z,E) 

d 



dE 



Ni{z,E), 
(B16) 
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